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ABSTRACT 

We show that a joint analysis involving independent cosmological probes at inter- 
mediate redshifts provides a remarkable cross-checking for the Hubble constant 
H which is competitive with the existing measurements of the local Universe. 
Although dependent on the physics at z ~ 1, this result is fully independent 
on the calibrations involved in the cosmic distance ladder. The cosmological 
probes to be considered here are: (i) Angular diameter distances of galaxy clus- 
ters through SZE/X-ray technique (0.14 < z < 0.89), (ii) the ages of old galaxies 
at intermediate redshifts (0.62 < z < 1.70), (iii) measurements of the Hubble 
parameter H(z) (0.1 < z < 1.8), and (iv) the baryon acoustic oscillation (BAO) 
signature (z = 0.35). By taking into account statistical plus systematic errors 
and assuming a flat ACDM cosmology (Hq and Qm as free parameters), our joint 
analysis provides Hq = 73.4+g} km s _1 Mpc" 1 (la) when only the first three 
probes are combined. By adding the BAO scale, we obtain Hq = 74.1+2 2 km 
s _1 Mpc -1 (la) which is a 3.0% determination of the Hubble constant at inter- 
mediate redshifts. Due to this special combination of tests, the present value of 
Hq is competitive with the latest determinations based on nearby Cepheids and 
SNe la [Riess et al. ApJ 730, 119 (2001)]. This value can be much improved 
in the near future when more and larger samples (with smaller statistical and 
systematic uncertainties) become available. The present result also suggests that 
the method proposed here can be useful to achieve the wished theoretical and 
observational convergence on the value of Hq. 

Subject headings: Cosmological test, Sunyaev-Zeldovich effect, Galaxy ages, cos- 
mic distances, Hubble constant 
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1. Introduction 

One of the most important observational quantities for cosmology is the Hubble con- 
stant, Ho, whose value determines the present day expansion rate of the Universe. Eight 
decades after the original Hubble measurement, the determination of H is still a very active 
subject (Tammann, Sandage & Reindl 2008; Mould & Sakai 2009; Freedman & Madore 2010; 
Tammann & Reindl 2011; Riess et al. 2009, 2011). Nowadays, several groups and missions 
are focusing their efforts on its precise determination since H works like a key to quantify 
many astronomical phenomena in a wide range of scales involving galaxies, clusters, and 
superclusters. It also plays an important role for several cosmological calculations as the 
physical distances to objects, the age, size, and the matter-energy content of the Universe 
(Peebles 1993; Peacock 1999; Freedman et. al. 2001; Jackson 2007). 

Currently, the more robust constraints on Ho are obtained from local tests (z <C 1). The 
method is based on the cosmic distance ladder through a combination of different calibrators. 
The basic measurements and strategies commonly adopt: Cepheids, tip of the red giant 
branch, maser galaxies, surface brightness fluctuations, the Tully-Fisher relation, and type la 
supernovae (for a review see Freedman & Madore 2010). Along the last decade, Hubble Space 
Telescope (HST) Key Project did a magnificent job for decreasing significantly the errors 
on Hubble constant (Friedman and Madore 2010). Recently, Riess et al. (2011) also used 
the HST observations to determine H from Cepheids and Supernovae (SNe la). Through a 
rigorous analysis of the statistical and systematic errors they obtained H = 73.8 ± 2.4 km 
s _1 Mpc -1 (la), corresponding to a 3.3% uncertainty. However, other authors also based 
on the HST observations have found H < 70 km s^ 1 Mpc -1 from a different treatment. 
Probably, the main difficult of the methods based on the cosmic distance ladder is the possible 
existence of Hubble bubble or other local effects that would affect such measurements (Jha 
et al. 2007; Conley et al. 2007; Neill et al. 2007; Sinclair et al. 2010). The possibility of 
an observational convergence in the near future has increased in the last few years, however, 
additional progress, say, for a determination with 1% uncertainty will require to rebuild some 
aspects of the cosmic distance ladder (Suyu et al. 2012). 

On the other hand, we know that the cosmologists desire accurate measurements of the 
H mainly to refine the constraints on the neutrino masses (^2m u ), density (fU), and the 
equation of state parameter (u) of dark energy based on CMB anisotropies data (Macri et 
al. 2006; Sekiguchi et al. 2010). Spergel et al. (2007) have shown that CMB studies can 
not supply strong constraints on the value of Ho on their own. This problem occurs due to 
the high degree of degeneracy on the parameter space (Tegmark et al. 2004) since different 
grouping of the parameters (Ho, Qm, ^a, etc.) produce the same prediction of the cosmic 
microwave background radiation (CMB) anisotropies. It may be circumvented only by using 
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a prior on H from independent measurements (Hu 2005). In this connection, Komatsu et 
al. (2011), used CMB, Supernovae (SNe la) and Baryon acoustic oscillations (BAO) data, 
to derive a model-dependent value of H . It should be recalled that the phenomenology 
underlying the CMB test works at very high redshift (z ~ 1070), during the decoupling 
between radiation and matter when the universe was only 370 million years old. Therefore, 
this determination of Hq involves a combination of physics and phenomena at very different 
scales and epochs of the cosmic evolution (low, intermediate and high redshifts), and, more 
important to the present article, due to the inclusion of SNe la data, the derived value of 
H is also somewhat dependent on the cosmic distance ladder. In such circumstances, it is 
particularly relevant to obtain accurate bounds on the value of H from physics relying on 
many different kinds of observations, and, as far as possible, to avoid the combination of 
phenomena occurring at very different epochs (like SNe la and CMB). 

In the last few years, some cosmological tests at intermediate redshifts (z ~ 1) have 
emerged as a promising technique to estimate H (Simon et al. 2005; Deepak & Dev 2006; 
Cunha, Marassi & Lima 2007; Lima, Jesus & Cunha 2009; Stern et al. 2010; Busti et 
al. 2012). The main advantage of these methods is their independence of local calibrators 
(Carlstrom et al. 2002; Jones et al. 2005). In principle, such methods would provide a cross- 
checking for local direct estimates which are free from Hubble bubbles, since the majority of 
galaxy clusters are well inside the Hubble flow. However, all these methods are dependent on 
the assumptions about the astrophysical medium properties, as well as, of the cosmological 
model adopted in their analysis. It has also been argued that such probes, individually, are 
not yet competitive with the traditional methods based on the cosmic distance ladder (see, 
for instance, Freedman et al. 2001). 

In this letter, we identify a remarkable complementarity involving four different cos- 
mological probes at intermediate redshifts (z ~ 1) whose combination provides a valuable 
cross-checking for the Hubble constant value. The first test is based on measurements of the 
angular diameter distance from galaxy clusters via Sunyaev-Zeldovich effect (SZE) combined 
with measurements of the X-ray flux (SZE/X-ray technique). It was suggested long ago by 
many authors, but only recently it has been applied for a fairly large number of clusters 
(Bonamente et al. 2006; Cunha et al. 2007; Holanda et al. 2012a; 2012b). The second one 
is the age estimates of objects like galaxies and quasars at intermediate redshifts. It became 
possible through the new optical and infrared techniques together the recent advent of large 
telescopes (Alcaniz & Lima 1999; Lima & Alcaniz 2000; Alcaniz et al. 2003; Deepak & 
Dev 2006; Lima et al. 2009). The third possibility arises from measurements of the Hubble 
parameter, H(z), from differential ages and baryon acoustic oscillations (Simon et al. 2005; 
Gaztanaga et al. 2009; Stern et al. 2010). As we shall see, the cooperative interaction among 
these three independent probes already reduce greatly the errors on the Hubble constant, 
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and, more interesting, located just on a redshift zone distinct both from CMB anisotropies 
and the methods defined by the cosmic distance ladder (Cepheid, SNe la, etc.). In principle, 
it is possible to derive tighter constraints on the possible values of H once the cosmology 
is fixed. However, in our analysis we marginalize over the cosmology by adding the fourth 
test which is represented by the BAO signature. The advantage of a combination with BAO 
(a scale independent of H ) was suggested earlier by Cunha, Marassi & Lima (2007), how- 
ever, restricted only to the Hubble constant determination based on the SZE technique. In 
what follows, we show that the distinct probes adopted here act in concert to predict tighter 
constraints on the value of H Q . 



2. Basic Equations, Probes and Samples 

Let us now consider that the Universe is described by a flat Friedmann-Robertson- 
Walker (FRW) geometry (in our units c = 1) 

ds 2 = dt 2 - a 2 (t) [dr 2 + r 2 (d6 2 + stn 2 6d(f))] , (1) 

whose evolution is driven by cold dark matter plus a cosmological constant (ACDM). 

In this background, the angular diameter distance (ADD), Da, can be written as (Lima 
& Alcaniz 2002; Lima et al. 2003; Holanda et al. 2010) 

_ , , n N 3000/i- 1 r dz' 

D , (2;A , nM) = _y o __Mpc, (2) 

where h = H /100 km s" 1 Mpc -1 (subscript denotes present day quantities), and the 
dimensionless function 'K{z!\ Vt M ) is given by H = [J1m(1 + z') 3 + (1 — ^m)] 1 ^ 2 - 

Bonamente and collaborators (2006) determined the ADD distance to 38 galaxy clusters 
in the redshift range 0.14 < z < 0.89 using X-ray data from Chandra and Sunyaev-Zeldovich 
effect data from the Owens Valley Radio Observatory and the Berkeley-Illinois-Maryland 
Association interferometric arrays. Assuming spherical symmetry, the cluster plasma and 
dark matter distributions were analyzed by using a hydrostatic equilibrium model account- 
ing for radial variations in density, temperature and abundances. The common statistical 
contributions for this galaxy clusters sample are: SZE point sources ±8%, X-ray background 
±2%, Galactic N H < ±1%, ±15% for cluster asphericity, ±8% kinetic SZ and for CMB 
anisotropy < ±2%. Estimates for systematic effects are as follow: SZ calibration ±8%, 
X-ray flux calibration ±5%, radio halos +3% and X-ray temperature calibration ±7.5%. 
Typical statistical errors amounts for nearly 20% in agreement with other works (Mason et 
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al. 2001; Reese et al. 2002, 2004), while for systematics we also find typical errors around 
+ 12.4% and - 12% (see table 3 in Bonamente et al. 2006). 

On the other hand, the age-redshift relation, t{z), for a flat (ACDM) model has also 
only two free parameters (H , Qm) 

t(z;h,Q M ) = g -i [™ _ dx (3) 

Note that for Qm = 1 the above expression reduces to the well known result for Einstein-de 
Sitter model (CDM, Q M = 1) for which t(z) = |i7 " 1 (l + z)~ 3 / 2 . As one may conclude from 
the above equation, limits on the cosmological parameters Qu and H (or equivalently h), 
can be derived by fixing t(z) from observations. Note also that the age parameter, T = H Q t z , 
depends only on the product of the two quantities Hq and t z , which are usually estimated 
from completely independent methods (Alcaniz & Lima 1999; Lima & Alcaniz 2000; Cunha 
k Santos 2004; Friaca et al. 2005). 

In this context, Ferreras et al. (2009) catalogued 228 red galaxies on the interval 0.4 < 
z < 1.3 using HST/ACS slitless grism spectra from the PEARS program thereby studying 
the stellar populations of morphologically selected early-type galaxies from GOODS North 
and South fields. The first subsample adopted here consists of six passively evolving old red 
galaxies selected from the original Ferreras et al. sample. The second data set of old objects 
is also a subsample of the Longhetti et al. (2007) sample which is composed by nine field 
galaxies spectroscopically classified as early-types at 1.2 < z < 1.7. They result from a near- 
IR spectroscopic follow-up of a complete sample of bright early red galaxies selected from the 
Munich Near-IR Cluster Survey (MUNICS) that provides optical (B, V, R, I) and near-IR (J 
and K') photometry. The galaxy ages were estimated by making use of stellar populations 
with the degeneracy between the age of the best-fitting models and the star formation time 
scale being broken through a mass- weighting of the masses. For both samples, the selected 
data set (eleven galaxies) provide the most accurate ages and the most restrictive galaxy 
ages. 

In Figure la, we display the sample constituted by the eleven data points chosen from 
two distinct subsamples of old objects. As discussed by Lima et al. (2009), we have added 
an incubation time with a conservative error bar for all galaxies. Such a quantity is an 
estimate of the amount of time interval from the beginning of structure formation process in 
the Universe until the formation time (tf) of the object itself (the time of the radiation era 
is negligible). It is also assumed that t inc varies slowly with the galaxy and redshift in our 
sample, and, in order to account for our ignorance about this kind of "nuisance" parameter 
(Fowler 1987; Sandage 1993), we have associated a reasonable uncertainty, <J tinc . In what 
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follows, we consider that t inc = 0.8 ± 0.4 Gyr. We also combine statistical and systematic 
errors (for more details, see Lima et al. 2009). 

In Figs, lb, (lc), we compare the age of these old objects at intermediate redshifts with 
the predictions of the ACDM models for different values of the free parameters (Q M and h). 
It is also worth noticing that the present status of systematic uncertainties in this context is 
still under debate. Jimenez et al. (2004) studied sources of systematic errors in deriving the 
age of a single stellar population and concluded that they are not larger than 10% — 15% 
per cent. Others authors consider systematic errors around 20% (Percival & Salaris 2009). 
In the present study we have adopted 15% for all data. 

The third observational probe comes from H(z) data, obtained from differential ages 
of galaxies and radial BAO (Simon et al. 2005; Gaztanaga et al. 2009). Some years ago, 
Jimenez and collaborators suggested an independent estimator for the Hubble parameter 
(differential ages of galaxies) and used it to constrain the equation of state of dark energy 
(Jimenez et al. 2002, 2003). Later on, the differential ages of passively-evolving galaxies were 
used to obtain H(z) in the range of 0.1 < z < 1.8 (Simon et al. 2005), and this sample 
was further enlarged by Stern et al. (2010). In addition, Gaztanaga et al. (2009) took the 
BAO scale as a standard ruler in the radial direction (Peak Method) thereby obtaining 
three more additional data: H(z = 0.24) = 79.7 ± 2.7, H(z = 0.34) = 83.9 ± 3.2, and 
H(z = 0.43) = 86.5 ± 3.5, which are model and scale independent. Now, by comparing the 
theoretical expression, H(z) = H [J1m(1 + z ) 3 + (1 — ^m)] 1 ^ 2 , with the observational data 
the corresponding bounds can be readily derived. We remark that the relative age difference 
(the key to the method) is only of the order of 2 — 3% (Stern et al. 2010). However, to be 
more conservative we are assuming here systematic errors of 8% for the relative age difference 
and radial BAO data. 

As we shall see next, a joint analysis based on the above 3 different probes already 
provides tight constraints on the value of Hq. However, as shown by Cunha et al. (2007), 
the analysis of T>a data (from SZE/X-ray technique) leads to more stringent constraints 
on the space parameter (Q M ,h) when combined with the BAO signature (Eisenstein et al. 
2005; Percival et al. 2010), and the same happens when the ages of high redshift objects are 
considered (Lima et al. 2009). Therefore, instead to fix a definite flat ACDM cosmology (for 
a given value of Om), it is natural to leave Qm free which will be more accurately fixed by 
adding the BAO signature as a fourth probe to the complete joint analysis performed here. 

It is widely known that the BAO peak (detected from a sample of 46748 luminous red 
galaxies selected from the SDSS Main Sample) is predicted to arise precisely at the measured 
scale of 100 /i" 1 Mpc. It is interpreted as a consequence of the baryon acoustic oscillations 
in the primordial baryon-photon plasma prior to recombination. Such a measurement can 
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be characterized by the dimensionless parameter 



A 



n 



1/2 
M 



,1/3 



2/3 



0.469 ±0.017, 



(4) 



where = 0.35 is the redshift at which the acoustic scale has been measured, and T(z*) is the 
dimensionless co moving distance to z*. Note that the quantity given by (@J is independent 
of the Hubble constant and, as such, the BAO signature alone constrains only the Qm 
parameter. This property is very characteristic of the BAO signature thereby differentiating 
it from many others classical cosmological tests. 



3. Complementarity for H 

To begin with, let us consider a joint analysis based only on the combination of the first 
three probes, namely: (i) SZE/X-ray distances, (ii) the age of the oldest intermediate redshift 
objects, and (iii) the measurements of the Hubble Parameter H(z). Further, a complete joint 
analysis including the BAO signature from the SDSS catalog it will be performed. For both 
analysis, we stress that a specific flat ACDM cosmology has not a priori been fixed. 

All the observational expression adopted here have only two free parameters (Ii,Qm)- 
In this way, we perform the x 2 statistics over the Qm — h plane. Combining the three tests 
discussed above the x 2 value reads: 

2/ I n _ \^ (PAzi] p) — T^Aobs,i) 2 \^ (t(zj] p) — Unc — t b s ,j) 2 \^ (H(z k ; p) — H obs ^ k ) 2 

X \ Z \P) — _2 _i_ _2 , _.2 + 2-~t rr 2 -L rr 2 _|_ rr 2-^1 rr 2 -L rr ' 

i a V Aob3il -T a stat + a syst j a t obsJ + °t lnc + a syst fe a H obStk + "syst 

(5) 

where the quantities with subindex u obs" are the observational quantities, <Jv Ao s is the 
uncertainty in the individual distance, a sta t is the contribution of the statistical errors, a tinc is 
the incubation time error, a sys t are the contribution of the systematic errors for each sample 
added in quadrature and the complete set of parameters is given by p = (h, Ojw). 

In Fig. 2a we show dimensionless Hubble constant (h), versus the the matter density 
parameter (Qm)- Our analysis combining statistical and systematic errors from galaxy ages 
are in black lines with 68.3%, 95.4% and 99.7% confidence levels (c.l.). In the same way, the 
green lines represent the confidence levels for the SZE/X-ray sample whereas the correspond- 
ing constraints from Hubble parameter H(z) sample are show in orange lines. As should be 
expected, the constraints of each sample on {h, Qm) plane are too weak when individually 
considered. However, due to the complementarity among them, our joint analysis for these 
three samples predicts H = 73.4 ±4.7 km s" 1 Mpc" 1 (6.4%) and Q M = 0.290i^ (la c.l.) 
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for two free parameters with xLin = 41.42. The reduced values are Xred = 0.702 (including 
systematics) and xled — 1 ( no systematics). 

As remarked earlier, in order to put tighter constraints on the parameters h and Vt M we 
must add the BAO signature which is independent of the Hubble parameter and, as such, 
can break the degeneracy on the mass density parameter (Cunha et al. 2007, Lima et al. 
2009). 

In Fig. 2b, we display the contours on the space parameter obtained through a joint anal- 
ysis involving the combination of all four probes. The dashed lines are cuts on (h, VLm) plane 
from BAO signature. The red, green and blue contours are constraints with 68.3%, 95.4% 
and 99.7% c.L, respectively. This complete analysis provides H = 74.II33 km.s _1 .Mpc _1 
(4.5% uncertainty) whereas the density parameter is VLm = 0.278lo' O 2 8 f° r two free parame- 
ters with a xLn = 41-53 ( X 2 red = 0.704). 

In Figure lc, we show the likelihood function for the h parameter in a flat ACDM 
universe. Both curves were obtained by marginalizing on the matter density parameter 
(JIm)- The shadow lines are cuts in the regions of 68.3% and 95.4% probability. For the 
solid black line the BAO signature was not considered. The constraints for the black line 
are h = 0.734 ±0.031 (0.064) with la (2a), respectively. For the red line the BAO signature 
has been included. The constraints are h = 0.741 ± 0.022 (corresponding to 3% error in 
h) and 0.045 (6.1%) with la and 2a, respectively. For all these analyses, the statistical and 
systematic errors were added. 

In table 1, we compare the results derived here with the latest determinations of the 
Hubble parameter. Note that competitive constraints were obtained using only probes of 
intermediate redshifts, and, perhaps, more interesting, the present constraints are consis- 
tent with recent results provided by independent methods at different epochs of the cosmic 
evolution. 



Table 1: Limits to h for different methods and epochs 

Reference Method 



h (la) 



Komatsu et al. 2009 ... 
Freedman & Madore 2010... 
Riess et al. 2011... 
This paper... 
This paper... 



CMB 
Cepheid Variables 
Cepheid+SNe+Maser 
SZE/X-ray+Age+H(z) 



0.705 ± 0.013 
0.730 ± 0.045 
0.738 ± 0.024 
0.734 ±0.31 



Epoch 



Z rsj 1100 

z ~ 
z ~ 
0.1 < z < 1.8 



SZE/X-ray+Age+H(z)+BAO 0.741 ± 0.022 0.1 < z < U 
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4. Conclusions 

Several ongoing and future experiments (HST Key Project, Si/oES, PLANCK, Spitzer, 
GAIA, JWST, and others) are dedicated to measure more precisely the present day value 
of the Hubble constant. Numerous fields of research in Astronomy will be benefit with the 
planned improvement for the 2% measurement of H likely to be achieved in the present 
decade. 

In this letter, we have demonstrated that a joint analysis involving four independent 
cosmological tests at intermediate redshifts provides a 3% determination of H including 
statistical plus systematic errors. Our approach was separated in two distinct parts. Firstly, 
we have implemented a joint analysis involving only three probes, namely: (i) angular di- 
ameter distance from galaxy clusters (SZE/X-ray technique), (ii) a sample of old objects at 
intermediate redshifts, and (iii) measurements of the Hubble parameter H(z) (at z ^ 0). 
This analysis has shown that such tests act in concert thereby predicting an unexpected 
tight value of H . Secondly, we have added the BAO signature which is independent of the 
Hq thereby obtaining a more accurate value of the Hubble constant. 

It should be stressed that our 3.0% determination of H became possible only when the 
BAO signature was implemented in our joint analysis. In the framework of a flat ACDM 
model we have obtained h = 0.739 ± 0.022 (including statistical plus systematic uncertain- 
ties). As shown in table I, this value is not only consistent, but has the same precision of 
a recent H determination using nearby Cepheids and Supernovae (Riess et al. 2011). It is 
advocated here that the present determination of the Hubble constant based only on probes 
at intermediate redshifts provides a competitive-crossing checking for any determination of 
H . We also believe that the method presented here has a potential at least comparable 
to other independent approaches. The main advantage of the present treatment is that it 
does not rely on extragalactic distance ladder being fully independent of any local calibrator. 
Naturally, its basic disadvantage rests on the large and different systematic uncertainties ap- 
pearing in each probe when separately applied. However, the remarkable complementarity 
among the four tests works in concert thereby reducing greatly the possible degeneracy and 
uncertainties appearing in the present determination of the Hubble constant. 
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Fig. 1. — Old galaxies data at intermediate redshifts as a cosmic probe, a) The age-redshift plane 
and the total sample of galaxies. Black points for z < 1.0 and z > 1.2 correspond, respectively, 
to the Ferreras et al. (2009) and Longhetti et al. (2007) samples. The black, red and blue 
bars represent the statistical, statistical+incubation and statistical+incubation+systematic errors, 
respectively. These data points correspond to the 11 galaxies of our selected subsample (see the 
main text), b) Age-redshift relation (the effect of J2m)- The age of the Universe for some selected 
values of the density parameter. All plots were drawn by choosing the h parameter as the best fit 
of our analysis (h = 0.74). For smaller values of VLm the ages of the Universe for a given redshift 
increase, thereby accommodating the oldest selected objects, c) Age-redshift relation (the h effect). 
Dotted curves are the predictions of the cosmic concordance model (Qm = 0.28, Qa = 0.7) and 
different values of h. We see that for values smaller than h = 0.6 or bigger than h = 0.9 the curves 
predicted by the models move away from the data. 
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Fig. 2. — The determination of Hq at intermediate redshifts. a) Hq and the complementarity of 
three different probes. Confidence regions (68.3%, 95.4% and 99.7%) in the (CIm, h) plane provided 
by the SZE/X-ray + galaxy age + H(z) data including statistical plus systematic errors. The 
best-fitting values are h = 0.734 and VLm = 0.290. b) The BAO razor. Contours in the £Im — h 
plane from a joint analysis involving the SZE/X-ray + galaxy age + H(z) + BAO data. The 
contours correspond to 1,2, 3a confidence levels. The best-fitting model converges to h = 0.741 
and £Im = 0.278. c) The final values of Hq. Likelihood functions for the h parameter in a flat 
ACDM Universe. As indicated in the figure, the dashed (red) curve correspond to a joint analyses 
involving from SZE/X-ray + galaxy age + Hubble parameter whereas the solid black curve includes 
the BAO signature. The horizontal lines are cuts in the regions of 68.3% and 95.4% probability. 
By including statistical plus systematic errors the joint analysis performed with the four different 
probes provides h = 0.741 ± 0.022 (1<t). As can be seen from Table 1, this constraint on Hq is 
in agreement with the latest independent determinations of the Hubble constant (Freedman and 
Madore 2010; Komatsu et al. 2009; Riess et al. 2011). 



